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We use a novel Monte Carlo method to study the Mott transition in an anisotropic triangular lat- 
tice. The real space approach, retaining extended spatial correlations, allows an accurate treatment 
of non trivial magnetic fluctuations in this frustrated structure. Choosing the degree of anisotropy 
to mimic the situation in the quasi-two dimensional organics, k— (BEDT-TTF)2Cu[N(CN)2]-X, we 
detect a wide pseudogap phase, with anomalous spectral and transport properties, between the 'un- 
gapped' metal and the 'hard gap' Mott insulator. The magnetic fluctuations also lead to pronounced 
momentum dependence of quasiparticle damping and pseudogap formation on the Fermi surface as 
the Mott transition is approached. Our predictions about the 'bad metal' state have a direct bearing 
on the organics where they can be tested via tunneling, angle resolved photoemission, and magnetic 
structure factor measurement. 



The Mott metal-insulator transition (MIT), and the 
proximity to a Mott insulator in doped systems, are cru- 
cial issues in correlated electron systems PHI] • The Mott 
transition on a bipartite lattice is now well understood, 
but the presence of triangular motifs in the structure 
brings in geometric frustration [5] . This promotes in- 
commensurate magnetic fluctuations whose nature, and 
impact on the MIT, remain outstanding problems. 

The organic salts provide a concrete testing ground for 
these effects 0H]. The k-(BEDT-TTF) 2 Cu[N(CN) 2 ]- 
X salts are quasi two dimensional (2D) materials where 
the BEDT-TTF dimers define a triangular lattice with 
anisotropic hopping [5]. The large lattice spacing, ~ llA, 
leads to a low bandwidth, enhancing electron correlation 
effects, while the triangular motif disfavours Neel order. 
The X^Cli-zBra; family shows a MIT as x drops below ~ 
0.75 |10j . The metallic state is very incoherent above T ~ 
50K: the resistivity [TT] is large, > lOOmfkm, the optical 
response has non Drude character [TJ1 [T3] , and NMR [HJ 
ITS] suggests the presence of a pseudogap (PG). How these 
properties arise in response to magnetic fluctuations, and 
the crucial low energy spectral features in the vicinity of 
the Mott transition, remain to be clarified. 

We use a completely new approach to the Mott tran- 
sition, using auxiliary fields, that emphasizes the role 
of spatial correlations near the MIT. Our principal re- 
sults, based on Monte Carlo (MC) on large lattices are 
the following, (i) The interaction ({/)-temperature (T) 
phase diagram that we establish has a striking corre- 
spondence with k-BEDT in terms of magnetic transi- 
tion and re-entrant insulator-metal crossovers, (ii) At 
intermediate temperature, in the magnetically disordered 
regime, we obtain a strongly non Drude optical response 
in the metal, and predict a pseudogap (PG) phase over 
a wide interaction and temperature window, (iii) The 
electronic spectral function A(k, uj) is anisotropic on the 
Fermi surface, with both the damping rate and PG for- 
mation showing a clear angular dependence arising from 
coupling to incommensurate magnetic fluctuations. 



To provide a quick background, there have been several 
studies of the single band Hubbard model on a triangular 
lattice [T6H28] to model organic physics. Dynamical mean 
field theory (DMFT) has been the method of choice [20]- 
[22], usually used in its cluster variant (C-DMFT) [23H2i] 
to handle short range spatial correlations. The results de- 
pend on the degree of frustration and the specific method 
but broadly suggest the following: (i) the ground state is 
a PM Fermi liquid at weak coupling, a 'spin liquid' PI at 
intermediate coupling, and an AFI at large coupling|16F 
HHj, (ii) the qualitative features in optics jjjj and trans- 
port [2DJ are recovered, (iii) there could be a re-entrant 
insulator-metal-insulator transition with increasing tem- 
perature for a certain window of frustration [241 125j . (iv) 
the low temperature SC state could emerge 29 -31 from 
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FIG. 1: U — T phase diagram of the Hubbard model at 
i'/t — 0.8. The phases are paramagnetic metal (PM), param- 
agnetic insulator (PI), antiferromagnetic metal (AFM) and 
antiferromagnetic insulator (AFI). The AFM and AFI are 
not simple Neel ordered. PG indicates a pseudogap phase, 
metallic or insulating. There is no genuine magnetic transi- 
tion in two dimensions so our T corr indicates the temperature 
where the magnetic correlation length becomes larger than 
the system size 24 x 24. The MIT line is determined from 
change in sign of the temperature derivative of resistivity, i. e, 
dp/dT = 0. 
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Hubbard physics, although there is no consensus [32 . 

Surprisingly, there seems to be limited effort on the 
nature of spatial fluctuations, which could be significant 
in this low dimensional frustrated system. To clarify this 
aspect we study the single band Hubbard model on the 
anisotropic triangular lattice: 

H = Y *y c 1<t c jV - V Y n i + u Y ™*T n 4 (!) 

(ij)a i i 

We use a square lattice geometry but with the following 
anisotropic hopping: — —t when Rj — Rj = ±ioo or 
±yao, where ao is the lattice spacing, and t%j — —t' when 
Ri — Hj = ±(i + y)ao- We will set * = 1 as the reference 
energy scale, t' = corresponds to the square lattice, 
and *' = * to the isotropic triangular lattice. We have 
studied the problem over the entire t'/t window [0,1], 
but focus on t'/t = 0.8 in this paper, fi controls the 
electron density, which we maintain at half- filling, n = 1. 
U > is the Hubbard repulsion. 

We use a Hubbard-Stratonovich (HS) transformation 
that introduces a vector field m^r) and a scalar field 
fa(r) at each site [33] IS] to decouple the interaction. 
We need two approximations to make progress, (i) We 
will treat the mi and fa as classical fields, i.e, neglect 
their time dependence, (ii) While we completely retain 
the thermal fluctuations in nij, we treat fa at the saddle 
point level, i.e, fa -> (fa) = (U/2)((n % )) = U/2 at half- 
filling. With this approximation the half-filled problem 
is mapped on to electrons coupled to the field nij (see 
Supplement). 

H eff = Y ^i^-p-gE mroH- m i ( 2 ) 

where ft = fi—U/2. We can write He// = H e i{m.i}+H c i, 
where H c i = (t//4) ^ i m? . For a given configuration 
{mi} one just needs to diagonalise H e i, but the {mj 
themselves have to be determined from the distribution 

Tr eet e-^'e-^' 
F{mi} ~ JVmTr^e-^e-^ (3) 

Equation (2) describes electron propagation in the m; 
background, while equation (3) describes how the m, 
emerge and are spatially correlated due to electron mo- 
tion. The neglect of dynamics in the m^ reduces the 
method to unrestricted Hartee-Fock (UHF) mean field 
theory at T = 0. However, the exact inclusion of classi- 
cal thermal fluctuations quickly improves the accuracy of 
the method with increasing temperature. We will discuss 
the limitations of the method further on. 

Due to the fermion trace, P{rrij} is not exactly cal- 
culable. To generate the equilibrium {nij} we use MC 
sampling 35 37] . Computing the energy cost of an at- 
tempted update requires diagonalising iJ ei . To access 
large sizes within limited time, we use a cluster algo- 
rithm 38J for estimating the update cost. Rather than 



diagonalise the full -ff e // for every attempted update, 
we calculate the energy cost of an update by diagonal- 
izing a small cluster (of size N c , say) around the refer- 
ence site. We have extensively benchmarked this clus- 
ter based MC method [38 . The MC was done for lattice 
of size N = 24 x 24, with clusters of size JV C = 8 x 8. 
We calculate the thermally averaged structure factor 
^(q) = ~W Sij ( m i ' mj)e lq '( Ri ~ Rj ) at each temperature. 
The onset of rapid growth in 5(q) at some q = Q, say, 
indicates a magnetic transition. The electronic proper- 
ties are calculated by diagonalising H e i on the full lattice 
for equilibrium {m^} configurations. 

Figfl] shows the U — T phase diagram at t'/t = 0.8. 
Our low temperature result is equivalent to UHF and 
leads to a transition from an uncorrelated paramagnetic 
metal to an incommensurate AF metal with wavevector 
Q = Qi - {0.85tt, 0.85tt} at U cl - 4.0t. At U c2 - 4.5* 
there is a transition to an AF 'Mott' insulator with Q 2 ~ 
{0.87T, 0.87r}. The magnitude njj = |rrij| is small in the 
AFM and grows as U /* increases in the Mott phase. The 
existence of the AF metal, and the nature of order in the 
intermediate U/t Mott phase, could be affected by the 
neglected quantum fluctuations of the m^. 

Finite temperature brings into play the low energy fluc- 
tuations of the Dij. The effective model has the 0(3) 
symmetry of the parent Hubbard model so it cannot sus- 
tain true long range order at finite T. However, our an- 
nealing results suggest that magnetic correlations grow 
rapidly below a temperature T corr , and weak interplanar 
coupling would stabilise in plane order below T corr . This 
scale increases from zero at U = U c i, reaches a peak at 




FIG. 2: Temperature dependence of the resistivity for dif- 
ferent U in the neighbourhood of U c . The normalising scale 
is po = hco/ire 2 (see text). This is ~ 380/if2cm for the organ- 
ics. The U /t values are indicated on the curves. The inset 
shows the experimental transport gap in the Ch-^Brx fam- 
ily (in Kelvins), and our estimated transport gap. We used 
a fit U(x)/t — 6 — 1.35a: — 0.4a; 2 (see text) to reproduce the 
transport gap[TT] estimated from the experiments. 
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U/t ~ 6.5, and falls beyond as the virtual kinetic energy 
gain reduces with increasing U. 

We classify the finite T phases as metal or insulator 
based on dp/dT, the temperature derivative of the resis- 
tivity. The dotted line indicating the MIT corresponds 
to the locus dp(T, U)/dT = 0. In addition to the mag- 
netic and transport classification we also show a window 
around the MIT line where the electronic density of states 
(DOS) has a pseudogap. To the right of this region the 
DOS has a 'hard gap' while to the left it is featureless. 
The MIT line shows re-entrant insulator-metal-insulator 
behavior with increasing T near U ~ U&. 

We can attempt a quick comparison of the phase di- 
agram with that in the k-BEDT family. The primary 
hopping is t ~ 65meV, and t' /t ~ 0.8 [H] (a recent ab 
initio estimate suggests t'/t < 0.5 for k-CY). Fitting 
the transport gap in k-CI (see Fig[l] inset) suggests that 
U/t ~ 6.0 at x = 0. From our results this would indicate 
that T corr /t 0.05, at x = 0, i.e, T corr ~ 35K, not too 
far from the NMR inferred T c ~ 30K. 

Figj2]shows the resistivity p(T) computed via the Kubo 
formula for varying U/t. We first compute the planar 
resistivity (which has the dimension of resistance) and 
then compute the effective three dimensional resistivity 
of 'decoupled' layers (see Supplement) by using the c-axis 
spacing, cq. In the C\i- x Ht x family it is observed that 
the transport gap can be fitted to A(x) « 800 — 1000a; 
Kelvinjllj. We match this to the U dependence of our 
calculated gap, A(U)/t, and infer U/t\ x=0 ~ 6. The MIT 
occurs at x c ~ 0.75 and for us at U/t ~ 4.5. Fitting a 
quadratic form to U(x)/t to capture the measured trans- 
port gap, we get U(x)/t « 6 - 1.35a; - 0.4x 2 . The U/t 
range in Fig [5] corresponds roughly to x — [0,1]. Since 
t = 65meV, T = OAt is approximately 300K. 

Our resistivity is in units of po = Using cq ~ 29 A, 
Po ~ 380/iilcm, yields p ~ 60po ~ 25mf2cm at T ~ 0.4i, 
while experimental value is > lOOmfkm. The difference 
could come from electron-phonon scattering absent in 
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FIG. 3: Optical conductivity at T jt = 0.1 and 0.2 for U vary- 
ing across U c . At these temperatures the ct(cj) is non Drude 
even in the 'weakly correlated' case U/t ~ 4.0. The finite 
frequency peak evolves into the Hubbard transition at large 
U/t. For a rough comparison to organics, T/t = 0.1 =80K, 
u/t = 5 = 2500cm" 1 , and a/a = 0.1 = 265fi" 1 cm- 1 . 




FIG. 4: Density of states at T/t = 0.1, 0.2 for U varying 
across U c . The dip in the DOS deepens with increasing T for 
U/t < 4.8. For larger U/t the system slowly gains spectral 
weight with increasing T. 

our model. Limelette et aZ[5D] presented DMFT based 
resistivity result that compares favourably with experi- 
ments, but, apparently, involves an arbitrary scale fac- 
tor. Our re-entrant window 6U=0At near U c , inferred 
from thermally driven I-M-I crossover, is equivalent to 
8x~{).2. This is consistent with 5x ~ 0.2 in the Cli-sBr^ 
family 11 . The C-DMFT estimates of the re-entrant 
window varies widely, from <5t7 ~0.3t[25 to ~1.2t[53]. 

Fig{3] shows the optical conductivity ct(uj) at T = O.lt 
and T = 0.2i as U/t varied across the Mott crossover. 
Our first observation is the distinctly non Drude nature 
of cr(u) in the metal, U/t ^ 4.4, with da(u>) /dui\ u ^o > 0. 
The low frequency hump in the bad metal evolves into the 
interband Hubbard peak in the Mott phase. The change 
in the lineshape with increasing T is more prominent in 
the metal, with the peak location moving outward, and 
is more modest deep in the insulator. 

In the organics the experimenters have carefully iso- 
lated the Mott-Hubbard features in the spectrum by re- 
moving phononic and intra-dimer effects |13j . Since we 
have already fixed our t, t', U we have no further fitting 
parameter for a(u}). The measured spectrum at x ~ x c 
and T ~ 50 - 90K has a peak around 1500 - 2000cm- 1 . 
Using U c /t ~ 4.5 and T/t = 0.1 we get u pea k/t ~ 3.0, 
which translates to ~ 1500cm -1 . The magnitude of 
our <j((jj) at ujpeak is ~ O.ldo ~ 265£!~ 1 cm~ 1 , since 
(To = 1/po ~ 2650f2 _1 cm~ 1 . This is remarkably close 
to the measured value, ~ 280i7~ 1 cm _1 (Ref. [13] Fig. 3). 

While the characteristic scales in a(ui) match well with 
experiments, the experimental spectrum has weaker de- 
pendence on temperature and composition x. This could 
arise from the subtraction process and the presence of 
other interactions in the real material. Our result differs 
from DMFT [TJ], and agrees with the experiments, in 
that we do not have any feature at to = U/2. We have 
verified the /-sum rule numerically. 

The crossover from the bad metal to the insulator in- 
volves a wide window with a pseudogap in the electronic 




4 




FIG. 5: Top: Momentum dependence of the low frequency spectral weight in the electronic spectral function A{k,uj) at 
T/t = 0.1. k x , k y range from [-it, it] in the panels. Note the systematically larger weight near k = [n/2, n/2] and [— 7r/2, — 7r/2] 
and smaller weight in the segments near [n,0] and [0, n]. U/t = 4.2, 4.4, 4.6, 4.8, 5.0, left to right. Bottom: Magnetic 
structure factor 5*(q) for the auxiliary fields mi for the same set of U jt. The q x , q y range from [0, 27r]. Note the very weak and 
diffuse structure at U/t = 4.2 and the much larger and differentiated structure at U/t = 5.0. 



DOS, N(u>). One may have guessed this from the de- 
pleting low frequency weight in Figj4] makes this 
feature explicit. We are not aware of tunneling stud- 
ies in the organics, but our results indicate a wide win- 
dow, U/t ~ [4,5.3], where there is a distinct pseudo- 
gap in the global DOS. This suggests that the entire 
x ~ [1.0,0.35] window in the organics should have a 
PG. For U/t < 4.6 the dip feature deepens with in- 
creasing T, we have dN(0)/dT < (compare panels (a) 
and (b), FigjiJ, while for U/t > 4.6 we have a weak 
dN(0)/dT > 0. The PG arises from the coupling of elec- 
trons to the fluctuating . A large at all sites would 
open a Mott gap, independent of any order among the 
moments. Weaker rrii, thermally generated in the metal 
near U c \ and with only short range correlations, manages 
to deplete low frequency weight without opening a gap. 
Since the typical size (mi) increases with T in the metal, 
we see the dip deepening at U < U c . 

While the size of the determine the overall deple- 
tion of DOS near u) = and the opening of the Mott 
gap, the angular correlations dictate the momentum de- 
pendence of the spin averaged electronic spectral function 
A(k, ui) (see Supplement). 

Within 'local self energy' picture, as in DMFT, A(\s.,uj) 
should have no k dependence on the Fermi surface (FS). 
In that case we should have k independent suppression 
of A(k, 0) with increasing U/t. 

Fig|5j top row, shows maps of A(k, 0) for k x ,k y = 
[— 7T, 7r]7 at T/t = 0.1, as increasing U/t transforms 
the bad metal to a Mott insulator. The first panel at 
U/t = 4.2 (roughly a Br sample) shows weak anisotropy 
on the nominal FS while Figj4j(a) suggests that a weak 
PG has already formed. At U/t = 4.4, next panel, the 
weak anisotropy is much amplified and the weight in the 



[0, 0] — > [it, 7r] direction is distinctly larger. The next 
three panels basically show insulating states but without 
a hard Mott gap. Overall, the 'destruction' of the FS 
seems to start near [tt, — n], the 'hot' region, and ends 
with the region near [tt/2, ir/2], the 'cold spot'. We show 
data on the full A(k, uj) in the Supplement that indicates 
that with increasing U a PG feature forms at the hot 
spot while the cold spot still has a quasiparticle peak. 

Second row in Fig [5] shows the S(q) of the auxiliary 
fields at T/t = 0.1 for the same U/t as in the upper row. 
While there is no magnetic order we can see the growth 
of correlations centered around Q ss [0.857T, 0.857r] as U/t 
increases. The dominant electron scattering would be 
from k to k + Q, and the impact would be greatest in re- 
gions of the FS in the proximity of minima in |Vek|- The 
location of the hot spots on the FS, and the momentum 
connecting them, indeed correspond to this scenario. 

While we have a method that captures non trivial spa- 
tial correlations and its impact on electronic properties, 
we need to be cautious about some shortcomings, (i) Ne- 
glecting the dynamics of the misses correlation effects 
in the ground state of the metal and underestimates U c /t. 
(ii) It also misses the 'Fermi liquid' physics in the low T 
metal, but should be reliable in the T/t > 0.1 regime 
that we have focused on. (iii) There is potentially a 'spin 
liquid' insulator [37J HI] at intermediate U/t and t' /t = 1. 
We do not know of such results at t' /t = 0.8, but would 
prefer to emphasize our finite T results rather than the 
nature of the ground state. 

Conclusion: We introduced and explored in detail a 
method which retains the spatial correlations that are 
crucial near the Mott transition on a frustrated lattice. 
Using electronic parameters that describe the k-BEDT 
based organics we obtain a magnetic T c , metal- insulator 
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phase diagram, and optical response that reproduces the 
key experimental scales. We uncover a wide pscudogap 
regime near the MIT, and predict distinct signatures of 
the incommensurate magnetic fluctuations in the angle 
resolved photocmission spectrum. 
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SUPPLEMENTARY INFORMATION: 



Derivation of the effective Hamiltonian 



Our starting point is the Hubbard model 



T n 4 



We implement a rotation invariant decoupling of the 
Hubbard interaction as follows. First, one can write 

2 

where rij = + nq. is the charge density, Sj = 

I Ea /3 c ia&<xP c iP = 2oj is the local electron spin oper- 
ator, and riii is an arbitrary unit vector. 

The partition function of the Hubbard model is 



Z = J L»[c,c]e" 



S = / dr£(r) 

£ = ^2c la (T)d T c la (T) + H(t) 

ia 

We can introduce two space-time varying auxiliary fields 
for a Hubbard-Stratonovich transformation: (i) <j>i(r) 
coupling to charge density, and (ii) Ai(r)mj(r) = mj(r) 
coupling to electron spin density (Aj is real positive). 
This allows us to define a SU(2) invariant HS transfor- 
mation (see ref. [HIS]), 



(4? m 2 \ 



The partition function now becomes: 

_ /" t-t dcjdcjdfadmi (_ j 13 c(t)) 

£(r) = ^c iCT (r)<9 r c iCT (r) + i? (T) + £„ lt (0 l (r),m l (T)) 



nx 



^ + i^r^ + -j-j- - 2m, • si 



As discussed in the text, to make progress we need two 
approximations: (i) neglect the time (r) dependence of 
the HS fields, (ii) replace the field <j>i by its saddle point 
value (U/2)(rii} = U/2, since the important low energy 
fluctuations arise from the m,. Substituting these, and 
simplifying the action, one gets the effective Hamiltonian 



eff 



mf 

IT 



where p, = [i — U/2. For convenience we redefine rrij — > 
T^rnj, so that the is dimensionless. This leads to the 
effective Hamiltonian used in the text: 



The partition function can be written as: 

Z = J VmiTr c ^e~ m *ff 

For a given configuration {m^} the problem is quadratic 
in the fermions, while the configurations themselves are 
obtained by a Monte Carlo as discussed in the text. 

Optical conductivity 

The conductivity of the two dimensional system is first 
calculated as follows (ref. [3]), using the Kubo formula: 

°SSM = ~ E ^^MJxW 2 ^ - (e, - e a )) 

Where, the current operator J x is 

Jx = [i(cj i0 .Cj + £ !(7 - he) + t'(c\ !<7 Ci + $ + ^a- - be) 

i,(T 

The d.c conductivity is the lu — > limit of the result 

2 

above. ao= E ^-, the scale for two dimensional conduc- 
tivity, has the dimension of conductance. n a = f(e a ) is 
the Fermi function, and e a and |a) are respectively the 
single particle eigenvalues and eigenstates of H e ff in a 
given background {m^}. The results we show in the text 
are averaged over equilibrium MC configurations. 

The experimental results are quoted as resistivity of a 
three dimensional material. If we assume that the planes 
are electronically decoupled, as we have done in the text, 
then the three dimensional resistivity p^D can be esti- 
mated from the resistance of a cube of size L 3 . If the 2D 
resistivity is p^D = 1/ °~2D , the resistance of a L 2 sheet is 
just p2D itself. A stacking of such sheets, with spacing 
Co in the third direction, implies that the resistance of 
the L 3 system would be Rzd = P2DC0/L. By definition 
this also equals p^dL/L 2 = p^jj/L. Equating the two, 
P3D = PidCq. 

Spectral function 

We extract the thermal and spin averaged spectral 
function A(k,w) as follows. First, the retarded Greens 
function 

G a (^t) = -i9(t){{c ka (t),cl rT (0)}) 
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FIG. 6: The spectral function A(k,u}) at two k points on the FS that correspond to the highest and lowest value of A(k,0). 
We highlight the anisotropy over a range of U /t values, as the system evolves from a moderately damped metal to a pseudogap 
phase, and three temperatures. 



which can be simplified to 



where {\a}} are the single particle eigenstates and e Q are 
eigenvalues in a given {mi} background. In frequency 
domain, this becomes 



G ff (k,c) = ^ 



| (kale 



We have averaged the spectrum over the four k neigh- 
bours of the nominal 'cold' and 'hot' points of our 24 x 24 
lattice. This averaging reduces the anisotropy, so the true 
anisotropy would be greater than what we show here. 
Also note that at T — O.li, where the A(k, 0) in the text 
is shown, the spectral function has no peak at uj — ei- 
ther in the cold or hot spot. The pseudogap feature is 
visible all over the FS even at U/t = 4. 



From this: A a (k,Lu) = — — ImG (T (k, uj) is simply 
J2 a \(k<j\a)\ 2 5(uj — e Q ). We average this over spin ori- 
entations, er, and over thermal configurations. The k 
dependent weight at ui = is shown in Fig. 5 in the 
text. The full spectral function at the 'cold spot' and 
'hot spot', where A(k, is maximum and minimum, are 
shown, respectively, in the top and bottom panels in the 
figure [6j 
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